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SUMMARY 

A semianalytical approach is developed for the sensitivity analysis of 
linear unsteady aerodynamic loads. The semianalytical approach is easier to 
implement than the analytical approach. It is also computationally less expen 
sive than the finite difference approach when used with panel methods which 
require a large number of panels. The semianalytical approach is applied to 
an isolated airfoil in a two-dimensional flow and rotating propfan blades in 
three-dimensional flow. Sensitivity coefficients with respect to non-shape- 
dependent variables are shown for some cases. It is expected that the semiana 
lytical approach will be useful in aeroelastic design procedures particularly 
when mistuning is present and that it is potentially useful for shape sensitiv 
ity analysis of linear unsteady aerodynamics. 


INTRODUCTION 

The computation of derivatives of response quantities with respect to 
design parameters, known as sensitivity analysis, plays an important role in 
developing reliable and efficient procedures for design optimization of practi 
cal aerospace structures. Integrating the structural and aerodynamic design 
processes for developing better aerospace structures in an automated manner is 
gaining interest (ref. 1). Aerodynamic sensitivity analysis is necessary in 
order to make possible an efficient inter-disciplinary approach to the optimi- 
zation of aerospace systems. However, the interest in sensitivity analysis 
has been mainly confined to structural applications (refs. 2 and 3). A recent 
paper (ref. 4), stressing the need for aerodynamic sensitivity analysis, con- 
tained only one reference on the subject and none dealing with unsteady 
aerodynamics. The present paper attempts to partially fill this void. Dwyer, 
Peterson and Brewer (refs. 5 and 6) applied sensitivity analysis to boundary 
layer flow equations to compute the effect of various parameters, though not 
in the context of design optimization. Recently, Haftka, et al . (ref. 7) pre- 
sented a sensitivity analysis of steady aerodynamics for multidisciplinary 
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optimization. The scarcity of published literature is perhaps due to the very 
complex nature of the aerodynamic response phenomena of practical aerospace 
structures. This is particularly true of rotary wing systems such as propfans. 
The prediction of steady and unsteady aerodynamic loads on these systems is 
still a subject of intense research. In addition, the coupling of computa- 
tional aerodynamic codes with structural analysis and optimization programs is 
not straightforward. 

The generalized unsteady aerodynamic loads on a structure are, in gen- 
eral, functions of the flow conditions, the structural geometry, and the struc- 
tural motion. The derivatives of the generalized aerodynamic forces, with 
respect to the variables representing the structural geometry and motion, are 
useful in aeroelastic structural design procedures. Their computation consti- 
tutes the aerodynamic shape sensitivity analysis. Shape sensitivities are 
immensely useful in automated design of aircraft wings and rotating blades. 

For example, shape sensitivities can be used to compute optimal sweep distribu- 
tion on propfan blades. 

These sensitivity derivatives can be calculated using a straightforward 
finite-difference approach. However, experience in other engineering disci- 
plines suggests that a semianalytical approach may be more computationally 
efficient. The objective of the present paper is to present a semianalytical 
approach for sensitivity analysis of subsonic unsteady aerodynamics used in 
flutter analysis, for the purpose of computing the derivatives of the general- 
ized unsteady aerodynamic forces. 

As a first step, only non-shape-dependent variables, representing the flow 
conditions and the structural motion, are considered. Such variables include 
free-stream Mach number, vibration frequency, mode shape and rotational speed. 
In addition to aeroelastic optimization, the derivatives with respect to flow 
conditions and structural motion can also be useful in computational schemes 
for aeroelastic analysis (refs. 8 to 10) and in computing the derivatives of 
the flutter Mach number and flutter frequency with respect to structural design 
variables (refs. 11 to 13). 

In the following, the problem of numerical sensitivity analysis is intro- 
duced and the three possible approaches to the problem, (finite-difference, 
analytical, and the semianalytical) are described. The semianalytical approach 
is then described in further detail in the context of unsteady aerodynamic com- 
putations using panel methods. Finally, the application of sensitivity analy- 
sis to isolated airfoils in two-dimensional flow and to single rotation propfan 
blades in three-dimensional flow is presented along with some typical results. 


NOMENCLATURE 


A domain of integration 

AN abbreviation for the analytical approach 

bf parameter denoting coupling of first and second propfan normal modes 
C matrix of aerodynamic influence coefficients 
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FD abbreviation for the finite difference approach 

Hi (2) first order Hankel function of the second kind 
K kernel function 

k reduced frequency , (ub/U) 

[ normal to blade helical surface 

M free-stream Mach number 

0 matrix of generalized forces 

Q nm generalized force in n-th mode due to motion in m-th mode 
S propfan advance ratio (tip speed/axial speed) 

SA abbreviation for the semianalytical approach 

U free-stream fluid velocity 

U r flow velocity relative to a point on the blade surface 

v n complex amplitude of normal velocity fluctuation 

W normal velocity distribution on the lifting surface = 

4ir |C| v n e^ u ® / U in the case of propfan aerodynamics 

W m normal velocity distribution in m-th mode 

x a pitch-axis location 

a parameter 

(J 1-M 2 

A prefix denoting finite difference increment 
Ap pressure differential across the blade surface 

S 2 re 

Ap — 2 AP e 1uo called "pressure load parameter" 

v 

Ap m pressure differential across the blade surface due to motion in m-th 
mode 

S n normal displacement in n-th mode 

e nondimensional control point location 

9, 9q angular coordinates 
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£ dimensionless chordwise pressure station coordinate 
p Q air density 

x symbol for CPU times 

Q propfan rotational speed 

to blade vibration frequency 

to to/Q 

' prime denotes variation 


NUMERICAL SENSITIVITY ANALYSIS 

In the finite difference approach, the sensitivity coefficient f'(x) of 
a function f of a parameter x is evaluated by repeating the analysis after 
incrementing each parameter. That is, 


f ' (x) « 


f(x + Ax) - f ( x ) 
Ax 


( 1 ) 


where the subscript FD denotes the finite difference approach and Ax is 
the finite-difference step size. While only first order forward differences 
are considered, the comments apply equally to others. The finite difference 
approach is computationally expensive and accuracy is very much dependent on 
the selected step size Ax and round-off error accumulation. The proper 
selection of a finite difference step size is not trivial. However, the 
finite difference approach is easy to implement. 

The analytical approach to the problem of sensitivity analysis consists 
of evaluating the sensitivity coefficient by direct analytical differentia- 
tion. That is, 


f'Cx) -<P AN <«) <2) 

The subscript AN denotes the analytical approach. The analytical approach 
obviously results in the exact sensitivity coefficient. However, it is diffi- 
cult to implement because the function f(x) tends to be very complicated in 
practice and, in general, cannot be directly differentiated. Furthermore, many 
response functions of interest are computed by special purpose programs and 
analytical differentiation is impractical. In addition, in unsteady aerody- 
namic applications using panel methods, the cost of computing cp AN (x) is not 
much smaller than the cost of confuting f(x). Thus the cost of their imple- 
mentation is not offset by a substantial reduction in computational effort, as 
is often the case in structural applications. Hence, analytical methods of 
sensitivity analysis are not easily justified in unsteady aerodynamic 
appl i cations . 

The semianalyti cal approach consists of analytical differentiation of the 
original function with respect to an intermediate function, the derivative of 
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which is then evaluated by numerical differentiation. Thus, if f(x) = 
f(g(x)), then, 


f(x>»» (x) - ^ aiii.. + - Ax ) -9<x) 
r ^x; ~ cp SA ^x; - dg • Ax 


( 3 ) 


The subscript SA denotes the semi analyti cal approach. The semianalytical 
approach combines the efficiency of the analytical approach with the ease of 
implementation of the finite difference approach. The most beneficial combina- 
tion results if most of the computational effort of evaluating the function 
f(x) is contained in computing f(g) and most of the analytical complexity is 
contained in computing g(x). It can also be sometimes expected to result in 
more accurate derivatives than the finite difference approach. For example, 
if g(x) is nearly linear, then the SA approach results in much more accu- 
rate derivatives than the FD approach. 


THE SEMIANALYTICAL APPROACH TO THE SENSITIVITY 
ANALYSIS OF UNSTEADY AERODYNAMICS 

The integral equation expressing the relation between the upwash and the 
pressure distribution can be written as (ref. 14) 

W = [ Ap K dA (4) 

where W represents the normal velocity distribution on the lifting surface, 

Ap the distribution of the pressure differential between the upper and lower 
surfaces and K the kernel function. K may be formulated for steady, 
unsteady, two-dimensional, three-dimensional, incompressible, or compressible 
flow. Thus, the techniques are applicable to any of the above flow conditions, 
though the implementation presented herein is limited to two selected aerody- 
namic models. We consider unsteady compressible flow where structural motion 
is simple harmonic. The discussion, however, is also applicable to nonharmonic 
motion by replacing the real variable <o, representing frequency, by the com- 
plex variable s, representing the Laplace transform variable. The normal 
velocity W depends on the assumed motion of the lifting surface (structural 
motion) and is considered known. The kernel function for nonrotating systems 
depends on the flow Mach number and frequency of motion. For rotating systems, 
it also depends on the rotational speed. Given the coordinate position of the 
elemental area AA, the kernel function can be computed at any point on the 
lifting surface. Knowing W and K, the pressure distribution Ap is 
obtained by inverting the integral equation above. 

For aeroelastic applications, the primary interest lies not in the pres- 
sure distribution but in the generalized force contained in the matrix Q 
given by 


Q nm = f A P m 6„ d A (5 

nm Ja m n 

where S n is the normal displacement in the n-th mode shape and Ap m is the 
pressure distribution resulting from motion in the m-th mode shape. If steady 
aerodynamic displacements are neglected, the normal displacement 6 n is inde- 
pendent of the assumed vibration frequency, w and the flow conditions. 
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We seek the derivatives of the generalized aerodynamic forces, contained 
in the matrix 0, with respect to parameters representing flow conditions and 
structural vibratory motion but independent of the domain of integration, A. 
Because the variation of generalized forces with respect to the mode shape is 
simpler, mode shape sensitivity is considered separately. First consider 
parameters that do not depend on the mode shape. Let a be one such parame- 
ter. Differentiating equation (5) with respect to a. 


30 


nm 


3a 


3a 


S n dA 


( 6 ) 


Thus, the sensitivity of the generalized force is equal to the generalized 
force acting on the blade resulting from a pressure differential across the 
blade surface given by 3Ap m /3a. We refer to 3Ap m /3a as the "pseudo-pressure 
differential" corresponding to the parameter a. 


To obtain the pseudo-pressure differential, we differentiate equation (4) 
to get 

jp - f Ap dA = [ K dA (7) 

8a Ja 3a J/\ aa 

Equation (7) is an integral equation in the unknown 3Ap/3a and is identi- 
cal to equation (4) if the left hand side and 3Ap/3a in equation <7) are 
replaced by W and Ap respectively. Thus, the sensitivity of the pressure 
distribution is equal to the pressure distribution that gives rise to the nor- 
mal velocity distribution given by 


3W 

3a 




We refer to this velocity distribution as the "pseudo-upwash" distribution 
corresponding to the parameter a. The pseudo-upwash is analogous to the con- 
cept of the pseudo-load used in the sensitivity analysis of static structural 
response (ref. 3). 


To evaluate the pseudo-upwash, 3W/3a and 3K/3a must be computed. While 
it is often less difficult to evaluate 3W/3a analytically, the analytical 
evaluation of 3K/3a is generally time-consuming and difficult to implement. 
Therefore, the pseudo-upwash is evaluated by replacing the quantities 3W/3a 
and 3K/3a by their finite difference approximations, AW/Aa and AK/Aa 
respectively. Considering forward differences, 


and 


3W AW W(a + Aa) - W(a) 


3a 

~ Aa = 

Aa 

3k 

~ ^ _ 

K(a + Aa) - K(a) 

3a 

Aa - 

Aa 


( 8 ) 


(9) 


For sensitivity analysis with respect to mode shapes, we first note that 
the pressure distribution is a linear functional and the generalized force a 
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quadratic functional of the mode shape. Hence, we consider variations instead 
of derivatives. Denoting variations by primes, we obtain from equation <5), 


Q 


nm 


A Pm 5 n 

r m n 


dA + 



dA 


and from equation (4) , 


( 10 ) 


W 


m 



( 11 ) 


so that, once again, the variation Ap^ of the pressure distribution can be 

obtained by using W m as the pseudo-upwash distribution. Thus, for mode shape 
sensitivity analysis, we use equations (10) and (11) instead of equations (6) 
and (7) respectively. 


For simple harmonic motion of an airfoil, the upwash can be approximated 
as (ref. 14) 


dS 

V P) ■ U r ds 


m 


(P) + iw S (P) 
m 


( 12 ) 


so that 


w ; <p) = u r df (p) + iws ; (p) <i3) 

where U r is the flow velocity relative to the point P on the airfoil and s 
is the arc length along the blade surface at fixed span. 

A popular way to solve the integral equation given by equation (4) for 
the pressure distribution is the family of panel methods. In these methods, 
the singular pressure distribution is approximated by a panelwise constant- 
pressure distribution. The Kutta condition is implicitly satisfied by choosing 
specific control points. Panel methods can easily handle any planform and are 
most suitable for lifting surfaces of complex geometry. Fine paneling of the 
planform is necessary because of the rough approximation of panelwise constant- 
pressure distribution. 

Once a discretization scheme (such as a panel method) is devised, practi- 
cally all of the analytical complexity in computing the generalized forces Q nm 
is contained in the computation of the kernel function K. Thus, a semianalyt- 
ical sensitivity analysis scheme, using numerical differentiation for the 
kernel function derivative, is of considerable advantage in overcoming the 
implementation problems associated with the analytical approach. Also, the 
aerodynamic code used for analysis can be used for sensitivity analysis as 
well by simply replacing the real upwash distribution with the pseudo-upwash 
distribution and the real pressure differential with the pseudo-pressure dif- 
ferential. The approach adopted here is similar to the semianalytical approach 
popular in the sensitivity analysis of static structural response (refs. 15 and 
16) and implemented in some general purpose finite element programs (ref. 17). 
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The efficiency of the semianalytical approach in comparison to the finite 
difference approach depends on the size of the problem. The overall costs of 
computing the kernel function and of inverting the integral operator are 
roughly proportional to the square and the cube of the number of panels, 
respectively. Hence, if the number of panels is n, then the CPU time for the 
computation of the generalized forces can be expressed as 

t = ti n^ + x2^ (14) 

The first term represents the time required for inverting the integral 
operator and the second term the time required for kernel function evaluations. 
The semianalytical approach utilizes the same inverted integral operator 
(available from the solution of the corresponding analysis problem) for all 
derivatives, in contrast to the finite difference approach. Hence, the CPU 
times for the FD and the SA approaches of sensitivity analysis can be 
expressed as 


X FD = T i n ^ + x 2 n ^ (15) 

t SA = x 2 n2 (16> 

In pressure formulations, for a moderate number of panels, the second term 
of equation (14), representing the kernel function computation dominates the 
overall computational time, so that tp D and t sa are of the same order. 

However, as n becomes larger, the cost of the FD approach escalates more 
rapidly than that of the SA approach. Hence, when the number of panels is 
large, the efficiency gain of the SA approach over the FD approach can be 
substantial. A large number of panels is needed in the case of a complex geom- 
etry as that of a propfan, especially if mistuning is present (ref. 18). A 
further discussion of the computational cost of the solution of the integral 
equation is given by Clark and James (ref. 19). 


APPLICATIONS 

The semianalytical approach is applied to two unsteady aerodynamic models: 
(1) isolated airfoil in two-dimensional flow and (2) rotating propfan blades in 
three-dimensional flow. Sensitivity coefficients of generalized forces with 
respect to various parameters are calculated by both the FD and the SA 
approaches and compared. For the purpose of computing the finite difference 
derivatives, the step sizes were computed by a procedure described by Gill 
(ref. 20). This procedure seeks to minimize the total error (the sum of round- 
off and truncation errors) in the sensitivity coefficient f'(x) using the 
round-off error expected in the evaluation of function f(x). The round-off 
error in the evaluation of f(x) is estimated by the difference in the values 
of f(x) calculated in double precision and extended precision. 

In the following, the logarithmic sensitivity function is defined as x 
f'(x)/f(x). Thus, the sensitivity coefficient gives the absolute rate of 
change of the function f(x), whereas the logarithmic sensitivity function 
gives a measure of the relative rate of change of f(x). 
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ISOLATED AIRFOIL IN TWO-DIMENSIONAL 
COMPRESSIBLE FLOW 


The kernel functed for this unsteady aerodynamic model is the Possio ker- 
nel (ref. 21) given by 


K(x,5) - - gjp e lk(x 5> 


fk(x-C)/(3 


IX u (2) 


H 


<m|x|) kt 


(17) 


where M is the free-stream Mach number, k the reduced frequency (wb/U), « 
the blade oscillation frequency, b the blade semi-chord, U the free-stream 
flow velocity, and x and £ represent the dimensionless chordwise control 
and pressure stations respectively. The Possio kernel is valid for subsonic 
Mach numbers. Bland (ref. 21) gives a full description of the computational 
aspects of the Possio kernel evaluation. 


For the isolated airfoil, the parameters considered are the free-stream 
Mach number M, the reduced frequency k and the location of the pitch axis 
x a . The nominal point chosen for this configuration corresponds to M = 0.8, 
k = 0.10, and x a = -0.5 (i.e., the pitch axis is at the quarter-chord). The 
downwash and pressure are computed at 32 stations along the chord. The finite 
difference step sizes were selected by minimizing the estimated total error in 
the sensitivity coefficients of selected normal velocities, pressure differen- 
tials and generalized forces. The estimated error in the derivatives (ref. 20) 
was less than 1.5 percent in all cases. 


Figure 1 shows the variation of the generalized forces due to pitching 
motion and their sensitivity coefficients with respect to the parameter, k. 

Both the real parts (in phase with the motion) and the imaginary parts (in 
quadrature with the motion) are shown in figures 1(a) and (b) respectively. 

The axes on the right correspond to the function and those on the left to its 
sensitivity. The squares give the finite difference (FD) sensitivities and 
the crosses the semianalytical (SA) sensitivities. In all cases (except at 
zero reduced frequency), the FD and the SA coefficients were in very good 
agreement. It can also be seen that the variations of the generalized forces 
are consistent with the calculated values of the respective sensitivity coeffi- 
cients. Similar results were obtained for generalized forces due to plunging 
and for sensitivities with respect to M and x a . 

Figure 2 shows the logarithmic sensitivities, defined in the previous sec- 
tion, of the generalized forces with respect to k. At this nominal point, the 
imaginary part of the lift due to pitching is seen to be the most sensitive to 
reduced frequency and the real parts of the lift due to pitching as well as 
plunging the least sensitive. 


To compare the efficiency of the semianalytical method with the finite- 
difference method, the percentage saving of CPU time achieved by the semiana- 
lytical approach in figure 3 as a function of the number of chordwise stations. 
These CPU times were obtained in computing the sensitivity of generalized 
forces with respect to Mach number. The CPU time for the semianalytical 
approach ranges from 0.65 sec when the number of chordwise stations is 32 to 
63.7 sec when the number is 256. The CPU times were obtained on the Amdahl 
5860 computer using VS FORTRAN Version 2 compiler with no compiler 
optimization. 
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PROPFAN BLADES IN THREE-DIMENSIONAL 
COMPRESSIBLE FLOW 


The second unsteady aerodynamic model studied is based on a three- 
dimensional linear subsonic lifting surface theory in frequency domain and is 
developed by Williams and Hwang (ref. 22) for aeroelastic analysis of propfans 
(ref. 23). 


The aerodynamic sensitivity analysis described below is implemented in a 
computer program called ASTR0P3 at NASA Lewis Research Center. ASTR0P3 is an 
aeroelastic analysis program developed for aeroelastic analysis of single rota- 
tion propfans (ref. 23). An automated flutter search procedure (ref. 10), 
requiring the generalized force sensitivities, is also incorporated in ASTR0P3. 

The integral equation expressing the relationship between the upwash and 
pressure distributions on a thin propfan blade is given by 

H< p) - - | A 3p< p 0 > | [K< p , p 0 >] dA 0 (18) 

Note that the definition of the kernel function does not exactly corre- 
spond to that of the previous section. 

Williams and Hwang (ref. 22) discretized the above integral equation by 
splitting the blade into n rectangular panels within each of which Ap is 
assumed constant. One control point is assigned to each panel at the midspan 
and lOOe percent chord position. The best value for e was empirically found 
by Williams and Hwang (ref. 22) to be 0.85. The discretization results in the 
algebraic system of equations given by 

H = C Ap Or, Ap = C _1 W (19) 

where W isa vector of the values of W at chosen control points on each of 
the panels, Ap is a vector of the values of Ap on each of the panels, and C 
is a matrix of aerodynamic influence coefficients given by 



“ [K(P,,P n )l dA n 
8© L i 0 J 0 


( 20 ) 


where the subscripts i and j refer to the panel numbers. 


The computation of the pseudo-pressure differential requires the deriva- 
tive of the vector of the pressure load parameters, 3Ap/3a. 3Ap/3a is 
obtained by differentiating equation (19), 


3H _ r 3Ap 
3a 3a + 



or. 


3Ap r -l r 3W 

3a " U L 3a 



( 21 ) 
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For mode-shape-dependent parameters. 


( 22 ) 


9Ap r -l 8 H 
8 a " L 8 a 

because the aerodynamic influence coefficients are independent of mode shapes. 

For propfan blades, the parameters considered are the free-stream Mach 
number M, the blade vibration frequency co, the rotational speed Q, the loca- 
tion of the control point e and a parameter bf denoting the coupling of the 
first and second normal modes of the blade. The parameter bf defines the 
normal displacement of the blade mode shape as S(P) = bf Si ( P) + ( 1 -bf ) 62 < P ) 
where Si and 62 are the normal displacements due to motion in the first and 
the second blade normal modes respectively. Thus, the parameter bf repre- 
sents the coupling of the first two normal modes. In this formulation, the 
domain of integration A is dependent on the advance ratio. In order not to 
change the shape of the helical surface and thus maintain the non-shape- 
dependent nature of the variables, the advance ratio was artificially held con- 
stant while computing the derivatives with respect to the Mach number and the 
rotational speed. This limitation will be removed when the sensitivity analy- 
sis is extended to shape-dependent variables. The control point location e 
is also considered a parameter because its value is chosen empirically and the 
sensitivity of the generalized forces to e is of interest. 

The planform and the normal modes of the SR3C-X2 propfan blade (ref. 23), 
shown in figure 4, were used as an example. The normal modes analysis was per- 
formed using COSMIC NASTRAN considering centrifugal deformations but not steady 
aerodynamic deformations. The planform was divided into 72 panels. The finite 
difference step sizes were selected by minimizing the estimated total error in 
the sensitivity coefficients of selected normal velocities, influence coeffici- 
ents, pressure differentials and generalized forces. The estimated error in 
the derivatives (ref. 20) was less than 2.4 percent in all cases. 


Sensitivity To Frequency, Mach Number and Rotational Speed 

The computed sensitivities with respect to the vibration frequency of the 
generalized forces due to motion in the first mode are shown in figure 5 along 
with the generalized forces on the right vertical axis. Figures 6 and 7 simi- 
larly show the sensitivities with respect to Mach number and the rotational 
speed respectively. Once again, the FD and the SA sensitivities are in rea- 
sonable agreement. Both the FD and the SA sensitivities are sometimes 
affected by round-off error as evidenced by their wavy nature. This indicates 
that the step-size selection algorithm used (ref. 20 ) may be inadequate in some 
cases and an improved algorithm (ref. 24) should perhaps be employed. 


Sensitivity to Empirical Parameters 

An important application of sensitivity analysis is in judging the sensi- 
tivity of generalized forces to empirical parameters. Such sensitivity gives 
a measure of the confidence in the values selected for the empirical parame- 
ters. As an example of such an application, the sensitivities of the general- 
ized forces are calculated with respect to the location e of the control 
point on the panel at which the normal velocity due to structural motion is 
specified. Figure 8 shows the logarithmic sensitivity of generalized forces 
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with respect to the control point location. It is observed that the general- 
ized force could be expected to change by up to 3 percent for a 1 percent 
change in the control point location. This is considered satisfactory. 

Finally, for the 72 panels used here, the percent reduction in CPU time 
from the semi analyti cal approach over the finite difference approach was 8 to 
13 percent. Greater reductions are expected for larger numbers of panels. 


CONCLUDING REMARKS 

A semianalytical approach to sensitivity analysis of linear unsteady aero 
dynamics is presented. Applications to isolated airfoils and propfan blades 
have shown that the semianalytical approach can be implemented relatively 
easily and that the semianalytical approach does not suffer from any accuracy 
problems compared to the finite difference approach. Preliminary studies show 
that the semianalytical approach will result in substantial savings in computa 
tional time for sensitivity analysis when the number of panels in propfan 
aeroelastic analysis is large. This is typically the case when mistuning is 
considered. 

It is expected that the semianalytical approach will be useful for comput 
ing shape sensitivity derivatives of generalized unsteady aerodynamic forces. 
These derivatives are expected to be used in developing efficient analytical 
and design optimization procedures for complex aeroelastic structures. 
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FIGURE 1. - SENSITIVITY OF GENERALIZED FORCES 
TO FREQUENCY. 
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FIGURE 3. - EFFICIENCY OF THE SEMI -ANALYTICAL APPROACH 
OVER THE FINITE-DIFFERENCE APPROACH. 
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FIGURE 4. - THE SR3C-X2 PROPFAN MODEL IN THE LEWIS 8x6 FOOT WIND 
TUNNEL. 
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FIGURE 5. - SENSITIVITY OF GENERALIZED FORCES TO FREQUENCY. 
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FIGURE 6. - SENSITIVITY OF GENERALIZED FORCES TO MACH NUMBER. 
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